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Abstract 

An analogy is explored between a setup of three atomic traps coupled via tunneling 
and an internal atomic three- level system interacting with two laser fields. Within 
this scenario we describe a STIRAP like process which allows to move an atom 
between the ground states of two trapping potentials and analyze its robustness. 
This analogy is extended to other robust and coherent transport schemes and to 
systems of more than a single atom. Finally it is applied to manipulate external 
degrees of freedom of atomic wave packets propagating in waveguides. 



1 Introduction 



Exploring the wave nature of massive particles has become possible through 
the enormous experimental advances in the cooling of neutral atoms, ions, and 
molecules to temperatures where the de Broglie wavelength becomes compara- 
ble to or larger than optical wavelengths. These achievements have stimulated 
great interest into the field of quantum atom optics as an analogue of quan- 
tum optics with light [1,2]. A major objective within this field is to develop 
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elements for the manipulation of the spatial wavefunction of atoms or atomic 
ensembles, as beam splitters, mirrors, lenses, etc. Applications are broad, rang- 
ing from a fundamental interest in probing the wave nature of particles to the 
manipulation of neutral atoms for implementing quantum gates and to the 
construction of atom interferometers for precision measurements of physical 
constants or as inertial sensors. In all these crucial requirement is - 

as in quantum optics - to preserve the coherence of the matter wave. 



Of special interest to atom interferometers as well as to quantum informa- 
tion processing are concepts to trap or to guide atomic matter waves. For 
trapped atoms, the interaction with external lasers can be precisely controlled, 
a spreading of the wave packet can be inhibited in some or all spatial dimen- 
sions, and the effect of gravity can be compensated. Trapping and guiding of 
neutral atoms is usually based either on the interaction of the atom's perma- 
nent magnetic dipole moment with magnetic fields [3] or on the coupling of 
laser fields to the atom's induced dipole moment [4] . Arrangements of current- 
carrying wires [5,3], superpositions of standing light-waves with different fre- 
quencies (superlattices) [6], or appropriately shaped microlenses [7,8] allow 
to design and control a variety of potential shapes. Examples are Y-shaped 
guiding geometries to split a wave packet [9], cold atoms storage rings from 
guides forming a closed loop [3,5], or traps whose separation can be controlled 
in time [8,10,11]. As two traps are brought to a close distance and tunneling 
takes place, an atom initially located in one of them oscillates in a Rabi-type 
fashion between the two potentials. This is in close resemblance to a two-level 
atom interacting with a laser field, but in contrast the 'Rabi-frequency' is con- 
trolled via tuning the tunneling interaction. Such a process, if implemented 
correctly, is coherent as it does not introduce uncontrollable phases, and it in- 
deed allows for a simple realization of quantum bits and quantum gates [12]. 
This technique however is not robust under variations of the system param- 
eters and thus requires precise temporal control of the potentials. The same 
problem is of course present in optical two-level system, where for this rea- 
son a variety of robust techniques have been developed, which are based on 
controlling couplings in multi-level systems, and which are nowadays standard 
techniques in experiments. 



Here we will provide a theoretical analysis of atom optical analogues to three- 
level techniques, especially discussing processes reminiscent of stimulated Ra- 
man adiabatic passage (STIRAP, [13]) to coherently move atoms between 
traps and and coherent population trapping (CPT, [14]) to create spatial su- 
perpositions of atomic wavefunctions. We will furthermore provide simulations 
showing that this technique is not only applicable to trapped atoms, but also 
to wave packets propagating in guiding potentials. 
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2 Time-dependent trapping potentials 



To obtain an analogy between external degrees of freedom of an atom in a 
system of trapping potentials coupled via tunneling and an electronic three- 
level system coupled via the electric dipole-interaction with two laser fields, 
consider a linear arrangement of three atom traps. We assume strong confine- 
ment in the orthogonal directions, such that the dynamics can be restricted 
to the one- dimensional (ID) Hamiltonian 

^free = j dx ijj^x) 

where V{x, t) describes the potential consisting of three traps with tunable 
distance. At each time t, iffree can be diagonalized to obtain the instantaneous 
eigenstates. For large distance between the traps, these arc localized at the 
centers of the traps. In the general case, states i^iix, t), (f)M{x, t), and (I)r{x, t), 
localized around the centers of the left, middle, and right trap, respectively, 
can be constructed from suitable combinations of the eigenstates with lowest 
energy. Restricting to these states, we can expand ip{x) — J2a=L,M,R ba<f>{x) to 
obtain [15] 

Hiree^ -JLM{t)b{bM - JMR{t)b^MbR-^ Yl /ia(^)^l^a + CC. (2) 

^ a=L,M,R 

Here Jap{t) = — J dx(j)l^{x,t)HfT.cc{x,t)(f)f3{x,t) describes nearest-neighbor tun- 
neling and = — / dx(f)'^{x,t)HfT.cc{x,t)(j)a{x,t) are the on-site energies. 
Interactions of next-nearest neighbors have been neglected. Considering just 
a single atom and shifting the ground state energy, we arrive at the following 
Hamiltonian isomorphic to the Hamiltonian of a three-level system coupled 
by two laser fields in the rotating wave approximation [16,13]: 

H^-JLMm\L){M\ + |M)(L|) - JMR{mM){R\ + \R){M\) 

-Mt) - iiLmM){M\ - iixnit) - iiLmR){R\- (3) 

The couplings —Jlm and ~Jmr ai'e the analogies of the Rabi frequencies of 
the pump and the Stokes laser, and iim—I^l 7^ and Hr—^l correspond to 
the detuning from the single- and two-photon transition, respectively (compare 
Fig. 1). 

Optical three-level systems have been extensively analyzed. Exploiting the 
different possible configurations of detunings and variations of the Rabi fre- 
quencies gives rise to a large number of coherent manipulation schemes of the 
underlying three-level system, among them stimulated Raman adiabatic pas- 



ip{x)= / dxip\x)Hf,e^ip{x), (1) 
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Fig. 1. Illustration of the analogy between a system of three coupled trapping 
potentials arranged linearly and an atomic three-level system in A-configuration. 
The tunneling matrix elements J correspond to the optical Rabi frequencies, the 
detunings are given by the difference between the on-site energies. 

sage (STIRAP, [13]), coherent population trapping (CPT, [14]), and electro- 
magnetically induced transparency (EIT, [17,18]). The analogy to the system 
of three coupled traps, as demonstrated by the Hamiltonian Eq. (3), clearly 
suggests to explore the application of these effects to coherently manipulate 
the external degrees of freedom of a trapped neutral atom, given the ability 
to control the corresponding trap parameters. Such a control should be possi- 
ble in various trapping configurations as optical [7,8,19] and magnetic [3,5,11] 
microtraps as well as in optical lattices by exploiting superlattice techniques 
[6]. 



In the following we will especially refer to neutral atoms trapped in arrays 
of optical microtraps created by illuminating a set of microlenses with a red 
detuned laser beam [8], such that in each of the foci of the individual lenses 
neutral atoms can be stored by the dipole force. By illuminating the mi- 
crolenses by independent laser beams under different angles, it is possible to 
generate various sets of traps which can be approached or separated by chang- 
ing the angle between the two lasers [8] . Changing the angle allows to control 
the couplings between the traps, the on-site energy can be tuned by changing 
the laser intensities. The optical potential generated by a single laser passing 
through a single lens has a gaussian shape, i.e., 

Krap(a;) = -'l^oexp (-^"^^^a;^) = -"^oexp (^-^("^)^) ' (4) 

where Vq and uj^ are depth and frequency of the potential, respectively, m is 
the mass of the neutral atom and = ^Jh/mu^ is the ground state size. As 
several potentials of this type are superimposed, controlling the coherent and 
adiabatic evolution of a trapped atom is complicate as the potential depth 
doubles if two traps lie on top of each other. To avoid this problem, either 
techniques to suppress non-adiabatic excitations can be employed [20,12], or 
the light intensity can be reduced in an appropriate way as the traps are 
approached, e.g., by adding a blue detuned laser to produce a compensating 
extra potential [21]. Here we will focus on the later option and assume V{x) = 
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mm[T4rap(a: + aL{t)), Vtr^p{x), Vtri,p{x - aR{t))], where aL{t) and aR{t) fix the 
centers of the traps. For m^C(j^a|y^/(8Vo) <^ 1 the potential then consists of 
three concatenated harmonic traps. 

The Hamiltonian (3) neglects contributions from non-adiabatic couplings to 
excited vibrational states as well as direct couplings from the left to the right 
trap. In what follows we will take into account the full Hamiltonian (1) through 
a numerical integration of the ID Schrodinger equation to simulate the dy- 
namics of a neutral atom in the three-trap potential. 



2.1 STIRAP - robust shifting of atoms between traps 



For zero detunings, one of the eigenstates of Hamiltonian (3), the dark state, 
only involves the states localized in the left and the right trap: 

|D(e)) =cose|L) -sine|i?). (5) 

Here is the mixing angle which depends on the couplings through cos O = 
Jlm/ Jmr- To implement a robust method to move an atom from the leftmost 
to the rightmost trap using tunneling, the counter-intuitive STIRAP sequence 
can be applied: first the right and middle traps are approached and separated, 
and, with an appropriate delay time tdeiay, the same sequence is used for the 
left and middle trap [Fig. 2 (a)]. This changes the mixing angle from = 
to = 7r/2 [Fig. 2 (b)], and if the atom initially is located in the left trap 
and the process is adiabatic, then the state is at all times identical to the dark 
state. This moves the atom directly from \L) to [Fig. 2 (c)]. The spatial 
wavefunctions {x\D) of the dark state for various times during the approaching 
process are plotted in Fig. 2 (d). 

The advantage of such a STIRAP-like process, as compared to a direct trans- 
port via Rabi-type oscillations, is its robustness with respect to the variation 
of certain experimental parameters. As shown in Fig. 3 (a), the scheme works 
for a large range of delay times ^Deiay a-nd minimum distances dmin- A similar 
robustness is found for variations of, e.g., the duration t,. of the approach- 
ing/separation process and the time ti for which the traps are kept at con- 
stant distance, the only requirements being the adiabaticity of the process and 
the order of approaching and separating the traps. In an experimental real- 
ization, certainly a shaking of the centers of the trapping potentials provides 
an important source of decoherence. It can, e.g., be caused by a mechanical 
vibration of the microlenses, or by variations in the laser phases for optical 
lattices. Here we anticipate a periodic variation of the distance of the traps 
with frequency well below the trapping frequency, namely a;shake — ^^"^oJx- As 
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Fig. 2. (a) Approaching sequence for a STIRAP-likc process, (b) the evo- 
lution of the tunneUng Rabi frequencies —JiMit) and —JMR{t) between the 
left and the middle and the middle and the right trap, calculated from 
the tunneling splitting energy of two traps, together with the mixing angle 
6(t) = ajctan {JLM{t) / Jm nit)), and (c) the corresponding ground state popula- 
tions; The parameters are d^^,_^a = d^^a = 9 (maximal distances between traps), 
d^^a = d^j^a = 1.5 (minimal distances), t^^iVx = t^^LO^ = 300 (time used to 
approach/separate the traps), tf^'^uJx = tf^^Ux = (time at which the traps are 
at the minimal distance, and tfj^iay^x = 120 (delay between the the approaching 
processes), (d) The spatial wavefunction {x\D) of the dark state for various times. 
The gray line gives the potential at uixt = 0. 

Fig. 3 (b) shows, the transport efficiency is not significantly degraded even for 
shaking amplitudes on the order of a few percent of the minimal trap distance 
if the delay time is appropriately chosen. 



A parameter difficult to control in an experiment is the exact horizontal align- 
ment of the traps. If a slight tilt is present in the potential, then gravity 
will change the relative depth of the potential minima. In this case, as has 
been reported in Bose- Einstein condensation in a double-trap potential [22], 
after a sufficiently adiabatic evolution, the atom(s) will eventually be found 
in the trap with lower energy. To allow for a transport to the desired state, 
the evolution should be explicitly non-adiabatic [23]. Here wc will show that 
the STIRAP-like transport is within a large range not affected by gravity. To 
this aim we add a potential 

Aytiit(a:^) = ^huJxOiX (6) 

stemming from gravity to the dipolar trapping potential V{x). The parameter 
7 determines the slope of the ramp, and we will use 7 > 0, such that the right 
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Fig. 3. Robustness of the atom optics version of STIRAP, i.e., the transfer efficiency 
from \L) and \R), measured by the population pji = All parameters not 

varied in the figures are as in Fig. 2. In (a) the delay time tdeiay between the two 
approaches (horizontal axis) and the minimal distances between traps (vertical axis) 
are modified, (b) shows the transfer efficiency as a function of ioeiay (horizontal axis) 
and of the amplitude of a shaking ashake in the positions of the outer traps (vertical 
axis) with lo shake = W~'^ll>x- For ashakc > the shaking of the outer traps is in 
phase, for ashake < it is out of phase by tt. 
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Fig. 4. Transfer efficiency for an additional tilted potential \tiit{x) = 7 hux ax, 
parametrized by 7. (a) Temporal variation of the energy levels obtained from a 
diagonalization of Eq. (3) for the parameters as in Fig. 2 and 7 = 0.02. Black arrows 
correspond to diabatic crossings, gray arrows indicate points where the probability 
of non-adiabatic transitions is larger. In (b) the transfer efficiency is plotted as a 
function of the time tr needed to approach and separate the traps and of the tilt; in 
(c) the dependence of the transfer efficiency on the tilt of the potentials is compared 
for STIRAP for the parameters from (b) with uixtr = 300 and for the transfer via 
Rabi-type oscillations between two traps for UxU = 300 {uxU = 12 such that full 
population transfer occurs for 7 = 0) and LUxtr = 32 {uxU = 25) . 

trap is shifted up in energy with respect to the left one. For the parameters of 
our simulations, a value of 7 = 10~^ corresponds to a difference in the potential 

energy of 3 ■ 10~^ hujx between the outer traps at the minimal distance. For 
7^1 such a tilt affects only the on-site energies —fiais in the Hamiltonian 
(3). In the picture of an optical A-system this corresponds to a shift from 
the one- as well as from the two-photon resonance. In this case there exist no 
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adiabatic path from \L) to [24]. This is exemphfied in Fig. 4 (a), which 
shows the energies of the three eigenstates of the Hamihonian Eq. (3) for the 
parameters of the STIRAP sequence from Fig. 2 and 7 = 2 x 10^^. To obtain 
the transport from the left to the right trap, the process has to be designed as a 
combination of diabatic (black arrows) and adiabatic (grey arrows) processes. 
However, the conditions to obtain a diabatic crossing at the points indicated 
by the black arrows are usually fulfilled, such that the transfer efficiency is 
dominated by the adiabaticity requirement. For this reason, for a given 7 the 
fidelity improves if the approaching and separating of the the traps is made 
slower as can be seen from Fig. 2. As should be stressed again, this is in 
contrast to the Rabi-type transport between two traps, where a faster process 
has larger fidelity of the atom ending up in the initially empty trap. Fig. 3 
(d) shows a comparison of the two schemes: next to the transfer efficiency for 
STIRAP for LOxtr = 300 as a function of the tile 7 two curves for Rabi-type 
oscillations are shown for a slow {uJxtr — 300) and for a fast {uJxtr — 32) 
approaching process. 



2.2 CPT-like and EIT-Uke effects 



From the isomorphism between the Hamiltonians, it is obvious that also other 
processes from three-level optics can be exploited here. The approach sequence 
can be modified to create spatial superposition states with maximum atomic 
coherence by evolving the mixing angle from = to © = 7r/4, correspond- 
ing to a delayed approach of the left trap, but to symmetric separation. This 
process, reminiscent of coherent population trapping (CPT), is similarly ro- 
bust to the variation of parameters as the STIRAP process, as long as the 
symmetry of the separation process is maintained. A detailed analysis can be 
found in [25], where also an EIT-like process is described. 



2.3 Effects of atom-atom interaction 



To arrive at the Hamiltonian of Eq. (3), the system has been simplified through 
a restriction to the lowest energy levels and to only a single atom. Our nu- 
merical simulations, which took into account excited states, show that for 
adiabatic processes and the potentials considered here the former simplifica- 
tion is justified and allows to reproduce effects known from A systems. Excited 
states however can be exploited explicitly by starting with the atom in a state 
different from the ground state [25]. On the other hand, also the restriction 
to a single particle can be released. For sufficiently low temperatures, interac- 
tion between bosonic atoms is dominated by s-wave scattering, and restricting 



8 



„ -10 

est 

1 " 
5 10 





(b) 




(e) 


1.0 - 








• 


t-0.8 - 


• 


4 


• 

• 


• 


CD 

^0.4- 


• 


















0.0 = 



-10 



10 -10 



10 -IS 
ct X (partiefe 1 ) 



10 -10 



10 




200 400 600 800 
Tin)et<t^ 



Fig. 5. Robust and coherent transport of a hole from the left trap to the right 

trap in a system of three traps filled with two bosonic ^^Rb atoms with scattering 
length at = lOGao- Shown are plots of the two particle probabilities |V'(a;i,a;2)P 
at four different times (a-d) as indicated by the arrows in (e). The initial state 
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again to states ^q, the Hamiltonian describing the system is modified as fol- 
lows [26,27]: 

n^HFree+l ^ UXb%k- (7) 

^ a=L,M,R 

Here Ua — 4:Trhdsc f dx\(})a{x)\^, where dgc is the ID scattering length which 
can be changed via changing the orthogonal confinement or exploiting a Fesh- 
bach resonance. If \Uq\ is sufficiently large to separate sectors in energy space 
with different particle number within one trap, then interaction allows to move 
coherently more than one atom at once, or to generate 'Schrodinger cat '-like 
states through a CPT sequence. On the other hand, starting from a system 
of three traps and two bosonic atoms initially in different traps allows to co- 
herently and robustly transport the 'hole', i.e., the empty site. An optical 
analogue of such a system has been studied in [28] , where coherent population 
trapping has been analyzed for two electrons with aligned spins in a three-level 
system. Also in this case a dark state exists which can be interpreted as the 
dark state of a 'hole'. A similar effect can be achieved in the atom optical sys- 
tem. As an example. Fig. 5 demonstrates the corresponding STIRAP process 
which moves the hole between the outer traps. 



3 Manipulation of matter waves in guiding structures 



In the previous part we have, in close analogy to the three-level processes for in- 
ternal atomic states, manipulated the external wavefunction of a trapped atom 
by a temporal variation of the coupling between traps. We will demonstrate 
that similar methods allow to manipulate an atomic wave packet propagating 
in an appropriately designed fixed guiding structure. We will assume a system 
of three waveguides oriented in the direction of the y axes, with |/-dependent 
distances (see Fig. 6 (a) for an example), and a corresponding Hamiltonian 



9 



-f^free = {Px +P^)/2m + V{x,y). Now, instead of considering the cigcnstates 
of the ID potential for each fixed time t. we can compute eigenstates at each 
position y, and as before combine the states with lowest energy to states 
4>a{x,y), « £ {L,M,R} localized around the center of each waveguide. The 
full wavefunction can then be decomposed as ip{x,y,t) — J2a(^a{y,t)<Pa{x,y), 
and inserting this expression into the Schrodinger equation gives the following 
equation for the evolution of the coefficients Ca{y,t) [9]: 



Here Hadiy) = J dx (f)*j{x,y){p'^/2m + V{x,y))(f)j3{x,y) are the Hamiltonian 
matrix elements for fixed position y and Kapi^y) = — / dx (j)*^{x,y)dy(()p{x,y) 
and Papiy) = — J dx 4>a{x, y)dy4>i3{x, y) are the kinetic and potential couplings, 
respectively. 



For time-dependent trapping potentials, the STIRAP or CPT sequence was 
induced by the counterintuitive temporal ordering of the approaching and 
separation processes of the traps. Similarly, in the case of waveguides we will 
apply such sequences in space. Thus, to obtain a STIRAP-like transport with 
the atom initially located in the left arm, first the right guide is approached to 
the middle one, then, with an appropriate delay, the tunneling is also switched 
on between the middle and the left guide. Finally, tunneling is turned off in 
the same order: first between the right and the middle and then between the 
middle and the left guide. For a CPT-like process to split an atomic wave 
packet coherently between two waveguides, the approaching sequence has the 
same counterintuitive order, but the separation is symmetric, see Fig. 6 (a). 
Now the additional coupling terms in Eq. (8) make the evolution more com- 
plex. Expanding Cq, into plane waves with momentum hk leads to a diagonal 
/c^-proportional term which accounts for broadening of the wave packet and 
to a term proportional to k Ka/s which induces velocity-dependent couplings 
between the waveguides. Furthermore, a velocity-independent modification of 
the couplings is introduced through the potential couplings Pap- 



To take again into account further effects beyond this illustrative approxima- 
tion in order to evaluate the performance of such processes we have numerically 
integrated the full 2D Schrodinger equation. We assume to initially have an 
atomic wave packet located in the left arm, with a gaussian profile in the direc- 
tion of the waveguide corresponding to mean momentum (ky) and momentum 

spread Aky = k,. {k,. = ^2mLJr/h; for the simulations cjr = ^x/S). In the 
transverse direction the wave packet corresponds to the ground state of the 
potential. Fig. 6 (b-d) shows an example of the time evolution in a structure 
which generates a sphtting of the wave packet through a CPT-like configura- 
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Fig. 6. (a) Contour plot of the waveguide potential used to split the wave packet 
incident in the left arm into a superposition of packets traveling in the left and 
the right arm using a CPT-like process. The contour line shown corresponds to 
a distance of 3/2 ground state widths away from the centers of the waveguides. 
(b)-(d) show density plots of \ip{x,y)\'^ at times tOxt = 20, ujxt = 60, Uxt = 120 
for a wave packet with mean momentum (ky) = 3.5kr and initial width Aky = kr- 
The minimal distance between waveguides is admin = 1-5. (e) The relative atomic 
fractions leaving the setup through the upper and lower exits of the structure as 
the minimal distance between waveguides is modified, (f) As (e), but for different 
mean velocities {ky). 

tion. During the evolution, the wave packet strongly broadens in the direction 
of propagation, but is nevertheless nearly equally split between the left and 
right outgoing arms with a negligible amount of reflection. This splitting is 
again relatively robust with respect to the parameters describing the poten- 
tial, provided the symmetry of the splitting is maintained. Fig. 6 (d) shows 
the change of the atomic fractions in the exit ports of the setup as the minimal 
distance of the waveguides is varied. The process is not as perfect as its coun- 
terpart in traps due to the additional couplings present here. Especially the 
velocity-dependent couplings modifying the desired CPT-like process play an 
important role: the larger the mean velocity, the stronger the deviation from 
the equal splitting, as can be seen from Fig. 6 (e). 



4 Conclusions 



We have studied the manipulation of the external wavefunction of an atom in 
a dipole potential consisting of three traps whose coupling can be changed as 
a function of time. As we have demonstrated, such a system, restricted to the 
lowest eigenstates, constitutes an analogon to the extensively studied system 
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of three internal atomic states coupled via two external laser fields. This al- 
lows to apply concepts as STIRAP, CPT, or EIT to coherently and robustly 
manipulate the external atomic wavefunction. Such processes are of potential 
interest, e.g., to move around atomic quantum bits or to create superposition 
states for interferometry. In particular, we have analyzed the robustness of a 
STIRAP-like process allowing transport of the atom between trapping poten- 
tials, and we have also shown that coherent processes are possible for several 
interacting atoms. 

As a different setup, we have studied atomic wave packets propagating in 
waveguide potentials, where the time dependence of the trap distances is re- 
placed by a spatial variation of the distance between waveguides. Due to addi- 
tional, partially velocity-dependent couplings, the evolution is more involved 
and the transport or splitting processes are not as clean as in the case of traps. 
Still, a stronger robustness as for schemes relying only on Rabi-type tunneling 
between traps can be achieved, as we have exemplified through demonstrating 
the coherent splitting of a wave packet between two arms, a scheme interesting 
for, e.g, interferometry. 
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